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An experimental and numerical investigation of the noise produced by high -subsonic and 
supersonic three-stream jets was conducted. The exhaust system consisted of externally- 
mixed-convergent nozzles and an external plug. Bypass- and tertiary- to-core area ratios 
between 1.0 and 2.5, and 0.4 and 1.0, respectively, were studied. Axisymmetric and offset 
tertiary nozzles were investigated for heated and unheated conditions. For axisymmetric 
configurations, the addition of the third stream was found to reduce peak- and high-frequency 
acoustic levels in the peak- jet-noise direction, with greater reductions at the lower bypass-to- 
core area ratios. For the offset configurations, an offset duct was found to decrease acoustic 
levels on the thick side of the tertiary nozzle relative to those produced by the simulated two- 
stream jet with up to 8 dB mid-frequency noise reduction at large angles to the jet inlet axis. 
Noise reduction in the peak- jet-noise direction was greater for supersonic core speeds than for 
subsonic core speeds. The addition of a tertiary nozzle insert used to divert the third-stream 
jet to one side of the nozzle system provided no noise reduction. Noise predictions are 
presented for selected cases using a method based on an acoustic analogy with mean flow 
interaction effects accounted for using a Green’s function, computed in terms of its coupled 
azimuthal modes for the offset cases, and a source model previously used for round and 
rectangular jets. Comparisons of the prediction results with data show that the noise model 
predicts the observed increase in low-frequency noise with the introduction of a third, 
axisymmetric stream, but not the high-frequency reduction. For an offset third stream, the 
model predicts the observed trend of decreased sound levels on the thick side of the jet 
compared with the thin side , but the predicted azimuthal variations are much less than those 
seen in the data. Also, the shift of the spectral peak to lower frequencies with increasing polar 
angle is over-predicted. For an offset third stream with a heated core, it is shown that including 
the enthalpy-flux source terms in the acoustic analogy model improves predictions compared 
with those obtained using only the momentum- flux. 


I. Introduction 

Exhaust- system noise emission continues to be an issue for commercial and military aircraft and is, therefore, 
the subject of ongoing research efforts as to its understanding, prediction and reduction. National and international 
regulations are imposing significant limitations on the noise footprint of aircraft operations, which are driving the 
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development of noise-reduction concepts. Assessment of proposed noise-reduction concepts involves the combined 
use of experimental testing and theoretical and numerical analysis. Detailed experimental measurements are necessary 
to screen potential noise-reduction concepts and to provide a database for the construction and validation of reduced- 
order models which can then be used in parametric studies to optimize the most promising designs. These reduced- 
order models can be purely (or semi) empirical [1], or statistically based methods, such as those based on an acoustic 
analogy (Kha varan, Bridges and Georgiadis [2], Goldstein and Leib[3] ). 

Future engine architectures may provide a third exhaust stream that will be available for potential noise 
reduction technologies. A third jet stream allows for additional geometric and parametric variation of the nozzle 
operation, and for an offset of the third stream relative to the core and bypass streams. The introduction of asymmetry 
into the flow field of the jet provides the potential for re-directing noise away from certain observer locations. 

Papamoschou and Debiasi [4] studied the effects of offsetting the fan stream of dual-stream supersonic jets. 
Their results showed a decrease in Mach wave radiation on the thicker side of the jet (created by the offset) and they 
attributed this reduction to increased mixing and a reduction of the potential core length with the offset fan stream. 
Flow and noise measurements in eccentric dual- stream jets by Zaman [5] seemed to confirm these results. 

A significant amount of subsequent experimental work was carried out by Papamoschou and co-workers and 
researchers at NASA ([6], [7], [8], [9], [10], [11]) testing and refining various offset concepts for dual- stream jets. 
Concepts investigated included the use of s-ducts, guide vanes and fan- stream wedges to create the desired offset. 
Computational studies ([12], [13]) were also carried out, based on the Reynolds-averaged Navier Stokes equations, to 
guide the design and selection of parameters used in these experiments and to help understand the impact of the 
different offset concepts on the flow. A model, based on Lighthill’s acoustic analogy formulation [14] , for the noise 
radiated by asymmetric, dual- stream jets was presented by Papamoschou and Rostamimonjezi [15] . 

Some of the concepts proposed and initially studied in the above works were systematically assessed in the 
Offset Stream Technology tests at NASA Glenn [16]. This test also extended some of these concepts, many of which 
were initially proposed for supersonic jets, to subsonic flows. The investigation focused on moderate (5) and high (8) 
bypass ratio exhausts. 

Initial experiments by Henderson [17] have shown the potential for noise reduction using a third stream under 
certain flow conditions. Papamoushco, Johnson and Phong [18] have carried out experimental work to study the noise- 
reduction potential of coaxial and offset three- stream jets. 

In this paper, results from an experimental and numerical investigation of the noise produced by high- subsonic 
and supersonic three-stream jets are reported. The exhaust system consisted of externally-mixed-convergent nozzles 
and an external plug. Various bypass- and tertiary- to-core area ratios were studied, with axisymmetric and offset 
tertiary streams, for heated and unheated conditions. 

In Section II the experimental approach is described. Results from the experiments are given in Section III. The 
acoustic analogy-based noise prediction method is briefly described in Section IV and results using this method are 
compared with experimental data for select cases in Section V. Conclusions and a discussion of the results are given 
in Section VI. 


II. Experimental Approach 

The experiments were performed in the Aero -Acoustic Propulsion Laboratory (AAPL) at the NASA Glenn 
Research Center shown in Fig. 1. The AAPL is a 20 m radius geodesic dome treated with acoustic wedges. The 
AAPL contains the Nozzle Acoustic Test Rig (NATR), which produces a 1.35 m diameter simulated forward flight 
stream reaching Mach numbers of 0.35 and contains the High Flow Jet Exit Rig (HFJER), a three- stream jet engine 
simulator capable of replicating most commercial turbo -fan engine temperatures and pressures [19]. 

Acoustic measurements were made with the far- fie Id array shown in Fig. 1. The array contains 24 
microphones located on a 13.7 m constant radius arc covering polar angles between 45° and 160°, where angles greater 
than 90° are in the downstream direction relative to the nozzle exit. All data were corrected for atmospheric absorption 
[20] and wind tunnel shear layer effects [21]. Data were acquired using V” Bruel and Kjaer microphones without grid 
caps, pointed directly at the nozzle exit. Microphone sensitivity and frequency response have been applied to all 
measurements. Narrowband results are presented as power spectral density on a one-foot lossless arc. One-third 
octave spectra are also presented on a one-foot lossless arc. 

The axisymmetric experiments used the externally-mixed, externally-plugged, convergent-nozzle system 
shown in Fig. 2 with the range of tertiary-to- core- area ratios (A t /A c ) and bypass-to-core-area ratios (Ab/A c ) shown in 
Fig. 3. All nozzle-system configurations used a core-nozzle exit diameter and area of 13.2 cm and 69.7 cm 2 , 
respectively, and a common bypass nozzle. The bypass-to-core-area ratio was varied by using core nozzles with 
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Figure 1. A photograph of the Aero-Acoustic 
Propulsion Laboratory (AAPL) showing the 
Nozzle Acoustic Test Rig (NATR) and the 
High Flow Jet Exit Rig (HFJER). 
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Figure 3. The nozzle design space used in the 
experiments. 
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Figure 2. The axisynmmetric-nozzle system 
used in the three-stream experiments. 

slightly different external contours which resulted in differences 
in the inner diameter of the bypass nozzle at the bypass-nozzle 
trailing edge. The tertiary area ratio was varied through a set of 
tertiary nozzles with the range of exit areas shown in Fig. 3. 

Two different approaches were used to create 
asymmetry in the third-stream flow: (1) the 

introduction of an offset duct upstream of the tertiary 
nozzle and (2) a tertiary nozzle insert that blocked a 
circumferential region of the third- stream nozzle exit 
(see Fig. 4). The offset duct, which was combined with 
the A t /A c =1.0 tertiary nozzle, produced a 0. 156” offset 
of the tertiary nozzle centerline relative to the 
centerlines of the core and bypass nozzles. The 
tertiary- nozzle insert blocked 228° of the third- stream 
nozzle exit resulting in A t /A c = 0.6. 

The conditions used in the experiments are 
shown in Table 1. The nozzle pressure ratio, NPR, is 
the ratio of the stagnation pressure of the jet to the 
ambient pressure. The nozzle temperature ratio, NTR, 




is the ratio of the stagnation temperature of the 
jet to the ambient temperature. Subscripts c, b, 
and t refer to the core, bypass, and tertiary 
streams, respectively. For heated core- stream 
conditions, NTRb = NTR t = 1.25. For unheated 
core conditions, the temperatures of the bypass 
and tertiary streams were also unheated. Jet 
conditions with NPR t = 2.10 produce inverted 
velocity profiles for the two outer jet streams. 

The experiments were conducted at simulated 
forward flight Mach numbers (Mfj) of 0.0 and 
0.3. For = 0.3, a simulated dual- stream jet, 
designated by NPR t = 1.0, was achieved by 
setting the third- stream conditions equal to those 
of the simulated flight stream. 

Particle Image Velocimetry (PIV) studies 

included two-component measurements in a streamwise plane with the hght sheet oriented along the centerline of the 
jet and stereo PIV measurements with a cross-stream orientation of the hght sheet. The two -component measurements 


(a) (b) 

Figure 4. Photographs of the offset nozzle systems using 
(a) the offset duct and (b) the third -stream nozzle inserts. 
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provided vector maps for up to 14 exit core diameters downstream of 
the plug tip. Stereo PIV measurements acquired data for one half of 
the jet and axial locations up to nine exit-core diameters downstream 
of the plug tip. 

The stereo PIV system was configured to provide cross-stream 
measurements of the 3 -component velocity field from the test article. 
The entire PIV system was mounted on a large traverse system to 
facilitate performing velocity plane surveys of the flow field. The 
entire cross-stream flow field could not be captured with sufficient 
spatial resolution to meet the test requirements. Hence, only the bottom 
half of the flow field was acquired. The Stereo PIV system employed 
two high-resolution (4008x2672 pixels) cameras, mounted in 
landscape mode, equipped with 180 mm focal length lenses and 8 mm extension tubes to provide a 526x272 mm 
(WxH) field of view. The PIV system was positioned so that the top edge of the field of view was approximately 25 
mm above the nozzle centerline. The cameras were mounted downstream of the model exit plane at nominally ±45° 
from the nozzle centerline. Stereo PIV calibrations were performed using a single plane target translated to 9 axial 
positions over a ±2 mm range. A 4th-order polynomial was used in the calibration and a calibration verification 
operation was employed to insure the calibration overlapped the laser light sheet plane. The measurement plane was 
illuminated with a dual head 400 mJ/pulse Nd:YAG laser system. The laser beams were formed into 1 mm by 550 
mm light sheets using cylindrical and spherical lenses. Both cameras were connected to a single computer system via 
a CameraLinkTM PCI card and the 400 frame pair data sequences were acquired and streamed to disk at a rate of 2 
frame-pairs/camera/sec. 

In order to facilitate a large field of view and high spatial resolution in the two-component PIV measurements, 
a four camera, 2x2 configuration was used. The 4008x2672 pixel stereo PIV cameras were used with their 4008 -pixel 
axis oriented vertically (portrait mode). The cameras were equipped with 180 mm focal length lenses and positioned 
so that their fields of view overlapped by 2.54 cm. A PIV calibration target was used to calibrate and register all four 
cameras. The physical registration of the four cameras was used in the setup of the vector processing grids in the top- 
left, top-right and bottom- left and bottom-right camera images so that no interpolation was required in the merging of 
the left/right vector maps. The final merged camera vector map covered an area of 355x560 mm (WxH). All 4 cameras 
were connected to a single computer system via two CameraLink PCI cards and the 400 frame pair data sequences 
were acquired and streamed to disk at a rate of 2 frame -pairs/camera/sec. 

Four independent seeding systems were required in this study: core flow stream, bypass stream, 3 rd stream and 
ambient flow. The heated core and bypass streams were seeded with 0.5 pm diameter alumina powder. A dispersion 
of the alumina seed material in ethanol was prepared using a pH stabilization technique [22]. The alumina/ethanol 
was dispersed in the flow well upstream of the nozzle using an air-assisted atomizing nozzle. The pH stabilization 
technique provides highly dispersed, unagglomerated seed particles in the flow. The tertiary stream was also seeded 
using the pH stabilized aluminum oxide dispersion. The ambient free -jet flow was seeded using a propylene glycol 
liquid seed material. Several fog generators were setup in the inlet tunnel to the free -jet - allowing nearly 18 m of 
mixing before entering the PIV measurement planes. 

The PIV image data were processed using multi-pass correlation with 64x64 pixel subregions on 32 pixel 
centers, followed by 32x32 pixel subregions on 16 pixels centers. Subregion distortion processing was also used to 
process the PIV data [23]. Subregion distortion was used to correct for velocity gradients across the subregion and to 
minimize the “peak- locking” effect, which is the tendency for the estimated particle displacements to preferentially 
concentrate at integer values. In the subregion distortion technique, the local velocity gradients surrounding the 
current correlation subregion are used to distort the subregion before the cross-correlation processing operation. 
Distorting the subregion yields correlation subregions with uniform particle displacements, and hence, reduces any 
bias caused by the velocity gradients. Typically two additional passes after the multi-pass processing are used with 
subregion distortion applied to refine the correlation peak estimates. Due to the oblique viewing of the model in the 
stereo PIV configuration, the nozzle was recorded in both the left and right camera views. The image of the nozzle 
corrupts the background in the image - leading to a loss of correlation in regions where the model is brightly 
illuminated by the laser light sheet. The Symmetric Phase Only Filtering (SPOF) technique was also applied in the 
data processing to mitigate any effects from the model being in the background of the images near the exit plane [24] . 
The final cross-stream velocity vector maps had 2 mm spatial resolution. The final 2-D streamwise velocity vector 
maps had a spatial resolution of 1.5 mm. Sequences of 400 velocity vector maps were acquired at each measurement 
station and ensemble averaged to provide first and second order statistics over the entire measurement plane. 
Chauvenet’s criteria was used to eliminate any outliers in the ensemble averaging process [25]. 
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Table 1 Experimental Jet Conditions 





PSD (dB) 


III. Experimental Results 

A. Axisymmetric Nozzle-System Results 

Three of the four tertiary nozzles produced discrete tones 
for some of the operating conditions in Table 1 . The tones, which 
were not trailing edge tones [26], were found to be the result of 
nozzle separations on the interior of the tertiary nozzles near the 
nozzle trailing edges. Roughening the nozzle surfaces with strips 
of sandpaper eliminated the tones as shown by the spectra in Fig. 
5 (where data for 140° have been offset by 5 dB for clarity). The 
surface roughening may have promoted boundary layer 
transition. However, the sandpaper resulted in a reduction of the 
third-stream area due and an associated reduction in broadband 
levels. Since the production of tones did not appear to introduce 
local broadband elevations, acoustic data were acquired with the 
untreated nozzles and tones subtracted from the resulting spectra 
before computing one-third octave spectra. 

The one-third-octave spectra obtained for the 
axisymmetric-nozzles at Mg = 0.3, NPR C = NPRb = 1.8, and 
NTR C = 3.0 are shown in Fig. 6 for a range of NPR t . The data 
for 100° and 140° have been offset 5 dB and 10 dB, respectively. 
For NPR t = 2.1, the velocities of the two outer streams are 
inverted since the tertiary- stream velocity is greater than that of 
the bypass stream. The introduction of the tertiary stream has little impact on low- and mid-frequency acoustic 
radiation at small and broadside angles (60° and 100°) to the jet for Ab/Ac = 2.5 and for A t /A c = 0.6. At these same 
angles, the introduction of the tertiary stream increases low-frequency acoustic radiation for Ab/A c = A t /A c = 1.0, with 
increases in NPR t resulting in increases in noise radiation. Peak frequency shifts are also noted with the addition of 

the third stream for broadside angles and Ab/A c =A t /A c = 1.0. In the peak-jet-noise direction, reductions in acoustic 

levels at peak and mid frequencies are achieved with the introduction of a third stream for all bypass-to-core and 
tertiary-to-core area ratios with the largest peak- frequency reduction occurring for Ab/A c = 2.5 with A t /A c = 0.6 and 
the largest mid- frequency reductions for Ab/Ac = 1.0. Increases in high-frequency noise occur at all observation angles 
and for all axisymmetric nozzle configurations when the third stream is introduced at the inverted velocity condition. 
Note that reductions in noise with the introduction of the third stream are accompanied by increases in ideal thrust. 
For Ab/Ac = 1.0, A t /A c = 1.0, and NPR t = 1.3, the ideal thrust increases by roughly 14% over that for the same bypass- 
to-core area ratio with NPR t = 1.0. 

The one-third-octave noise reductions, where noise reduction is given by the difference in the acoustic level for 
the two -stream jet (NPR t = 1.0) and that of the three-stream jet, are shown in Fig. 7 for the peak-jet-noise direction 
(140°). Positive values indicate the acoustic levels of the three-stream jet are lower than those of the two stream jet. 
The largest noise reduction (roughly 3 dB), achieved at mid frequencies, occurs for the smallest bypass-to-core area 
ratio, Ab/Ac =1.0 [see Figs. 7 (c) and (d)]. For equivalent exit areas on the core, bypass, and tertiary streams, up to 
2.5 dB reduction is achieved at low injection pressures (NPR t = 1.3) with no increase in high-frequency noise. For 
Ab/Ac = 2.5, maximum noise reduction is limited to 2 dB although this reduction is achieved at the peak frequency 
jet-noise frequency (100 Hz) for A t /Ac = 0.6 and NPR t = 1.5. While the results in Fig. 7 indicate the maximum 
achievable noise reduction is relatively insensitive to tertiary-to-core area ratio, results obtained at A t /Ac = 0.4 show 
little noise reduction for all bypass-to-core area ratios within the design space shown in Fig. 3 indicating the noise 
reduction trends in Fig. 7 cannot be applied to A t /A c <0.6. 

The results obtained at Mg = 0.0, NPR C = NPRb = 1.8, and NTR C = 3.0 for the Ab/A c = 2.5 and A t /Ac =1.0 
nozzle system are shown in Fig. 8 for a range of NPR t . A comparison of the data in Fig. 6 (a) with that in Fig. 8 
indicates that, for the same three-stream area ratios and operating conditions, greater mid-frequency noise reduction 
in the peak-noise direction is achieved at Mg = 0.0 than at Mfj = 0.3. At small and broadside angles to the jet, the 
impact of the third stream on acoustic radiation is similar for Mg = 0.0 and 0.3. 



Figure 5. The narrowband spectra acquired at 
the unheated conditions in Table 1, NPRt = 1.5, 
and the indicated polar observation angles. 
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Figure 6. The one-third octave spectra acquired at NPR C = NPRb = 1.8, NTR C = 3.0 and Mfj = 0.3 for 
the indicated observation angles. The data are for (a) Ab/A c = 2.5 and At/A c = 1.0, (b) Ab/A c = 2.5 and 


At/Ac = 0.6, (c) Ab/Ac = 1.0 and At/A c = 1.0, and (d) Ab/A c = 1.0 and At/A c = 0.6. 
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Figure 7. The noise reduction for NPR C = NPRb = 1.8, NTR C = 3.0 and Mfj = 0.3 at an 
observation angle of 140°. The data are for (a) Ab/A c = 2.5 and At/A c = 1.0, (b) Ab/A c = 2.5 and 
At/Ac = 0.6, (c) Ab/Ac = 1.0 and At/A c = 1.0, and (d) Ab/A c = 1.0 and At/A c = 0.6. 


The one-third octave spectra obtained for Ab/A c = 
2.5, A t /A c = 1.0, and Mg = 0.3 are shown in Fig. 9. The 
data for NPR C =1.5 and NPRb =1.8 are given in Fig. 9 (a) 
and for NPR C = 1.8 and NPR b = 1.5 in Fig. 9 (b). A 
comparison of the data in Figs. 6 (a) and 9 (a) shows that 
greater peak- and mid-frequency noise reductions are 
achieved in the peak-jet-noise direction for a high core 
velocity (NPR C =1.8) than for a low core velocity (NPR C = 
1.5) with the introduction of the tertiary stream. For all 
observation angles, the increases in low-frequency 
acoustic radiation with increasing NPR t at low core 
velocity do not occur at the higher core velocity. 
Additionally, increasing NPR t produces larger increases in 
high-frequency noise for the low core velocity than the 
higher core velocity. A comparison of the results in Figs. 
6 (a) and 9 (b) shows the addition of the third stream 
produces similar noise reduction characteristics in the 
peak-jet-noise direction for the same core velocity 
although, at the same NPR t , a slightly greater peak- 
frequency reduction is achieved for the higher fan velocity 
(NPRb = 1.8) than the lower fan velocity (NPRb = 1.5) and 
slightly greater mid- frequency reductions are achieved for 
the lower fan velocity than the higher fan velocity. At 
broadside angles (100°), the introduction of the third stream has a greater impact on noise reduction at peak- and mid- 
frequency noise for NPRb =1.5 than for NPRb =1.8. At small observation angles (60°), the impact of the tertiary 
stream on the resulting acoustic radiation appears to be insensitive to bypass- stream conditions. 



Figure 8. The one-third octave spectra acquired 
at NPRc = NPRb = 1.8, NTR C = 3.0 and M fj = 0.0 
for the indicated observation angles. The data 
are for Ab/A c = 2.5 and At/A c = 1.0. 
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Figure 9. The one-third-octave spectra acquired at Mfj = 0.3 and the indicated observation angles 
for Ab/Ac = 2.5 and At/A c = 1.0. The data in (a) were acquired for NPR C = 1.5, NPRb = 1.8, and NTR C 
= 3.0 and in (b) for NPR C = 1.8, NPR b = 1.5, NTR C = 3.0. 

B. Insert Nozzle-System Results 

The results obtained with the tertiary-nozzle inserts for Ab/A c = 2.5 are shown in Fig. 10 for NPR C = NPRb = 
1.8, NTR C = 3.0, and Mfj = 0.3. The data in Fig. 10 (a) and (b) were acquired with the insert located (and centered) on 
the side of the jet closest to the microphone array and on the side of the jet opposite to the microphone array, 
respectively. For both insert orientations, the introduction of the third stream at NPR t >1.7 increased acoustic 
radiation at all observation angles. For NPR t = 1.3, the introduction of the tertiary stream with the insert-nozzle system 
had only a slight impact on acoustic radiation. Data acquired for other third-stream area ratios [using inserts with less 
circumferential blockage than that shown in Fig. 4 (b)] produced trends similar to those in Fig. 10. 




(a) (b) 

Figure 10. The one-third octave spectra acquired at NPR C = NPRb = 1.8, NTR C = 3.0 and Mfj = 0.3 for 
the indicated observation angles using the third -stream nozzle insert and Ab/A c = 2.5. The data in (a) 
and (b) were obtained with the insert located on the side of the jet closest to the microphone array and 
on the side of the jet opposite to the microphone array, respectively. 

C. Offset Duct Nozzle-System Results 

The results for the offset-duct nozzle configuration obtained at Mg = 0.3, NPR C = NPRb = 1.8, and NTR C = 3.0 
are shown in Fig. 11. In Figs. 11 (a) and (c), the thick side of the tertiary nozzle [see Fig. 4 (a)] is on the side of the 
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Figure 11. The one-third octave spectra acquired at NPR C = NPRb = 1.8, NTR C = 3.0 and Mfj = 0.3 for 
the indicated observation angles using the tertiary -nozzle offset duct. The data in (a) and (b) were 
obtained with Ab/A c = 2.5 and in (c) and (d) with Ab/A c = 1.0. For (a) and (c) the thick side of the jet is 
closest to the microphone array. In (b) and (d), the thin side of the jet is closest to the microphone array. 


jet closest to the microphone array. In Figs. 1 1 (b) and (d), the thin side of the tertiary nozzle is on the side of the jet 
closest to the microphone array. For Ab/Ac = 2.5, increasing NPR t reduces peak- and mid- frequency acoustic radiation 
in the peak-jet-noise direction and increases low- and high-frequency radiation at small and broadside angles on the 
thick side of the jet [see Fig. 1 1 (a)] . On the thin side of the jet [see Fig. 1 1 (b)] , the introduction of the third stream 
at NPR t <2.1 has little impact on acoustic radiation in the peak- jet-noise direction and a similar impact on acoustic 
radiation to that for the thick side of the jet at small and broadside angles to the jet. Mid- frequency noise reduction in 
the peak-jet-noise direction on the thick side of the jet is slightly greater with the introduction of the third stream for 
Ab/A c =1.0 [see Fig. 11 (c)] than for the same NPR t and A\JA C = 2.5. At small and broadside angles to the jet, 
increasing NPR t results in greater increases in low- and mid- frequency acoustic radiation on the thick side of the jet 
for Ab/Ac =1.0 than for Ab/Ac = 2.5. On the thin side of the jet, the introduction of the third stream increases acoustic 
radiation at all frequencies and all observation angles for A\JA C = 1.0. For all area ratios and observation angles, the 
addition of the third stream at the inverted velocity condition results in elevated high-frequency levels. 

The one-third-octave spectra for the offset-duct nozzle system and supersonic core conditions (see Table 1) are 
shown in Fig. 12 for Ab/Ac = 2.5. All data have been acquired with the thick side of the nozzle closest to the 
microphone array. The data for 100° and 140° have been offset by 5 dB and 10 dB, respectively, in Figs. 12 (a), (b), 
and (c) and by 7 dB and 15 dB, respectively, in Fig. 12 (d). The circled regions indicate the presence of broadband 
shock associated noise. In the peak-jet-noise direction, noise reduction resulting from the addition of the third stream 
decreases with increasing NPRb for the same core and tertiary conditions and increases with increasing NPR C for the 
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(a) (b) 




(c) (d) 

Figure 12. The one-third octave spectra on the thick side of the offset-duct nozzle system with Ab/A c = 2.5 
for NTRc = 3.0 and M fj = 0.3. The data are for (a) NPR C = 2.1 and NPR b = 1.8, (b) NPR C = 2.1 and NPR b = 
2.1, (c) NPRc = 2.3 and NPR b = 1.8, and (d) NPR C = 2.3 and NPR b = 2.3. 



Figure 13. The one-third octave spectra on the 
thick side of the offset-duct nozzle system acquired 
with A b /A c = 1.0, NPRc = 2.3, NPR b = 1.8, NTR C = 
3.0, and Mfj = 0.3. 


same bypass and tertiary conditions. The largest noise 
reductions in the peak-jet-noise direction for the 
supersonic core conditions [see Figs. 12 (a) and (c)] are 
greater than those for the same bypass condition and 
subsonic core flow [see Fig. 11 (a)]. The addition of the 
third stream at NPR t <2.1 has little impact on acoustic 
radiation at small and broadside angles for all supersonic 
conditions and NPRb = 1.8. For NPR b = 2.3 [see Fig. 12 
(d)], the addition of the third stream increases 
broadband-shock noise levels. 

The results obtained for the offset-duct nozzle 
system and A b /Ac =1.0 are shown in Fig. 13 for NPR C = 
2.3, NPR b = 1.8, NTR C = 3.0, and M fj = 0.3. The data 
have been acquired with the thick side of the nozzle 
closest to the microphone array. A comparison of Figs. 
12 (c) and 13 show that greater noise reduction is 
achieved in the peak-jet-noise direction with the addition 
of the third stream for A b /A c =1.0 than for A b /A c = 2.5. 
At small and broadside observation angles, the addition 
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of the third stream increases broadband- shock noise levels for Ab/Ac = 1.0. 

The one-third-octave noise reductions in the peak-jet-noise direction on the thick side of the offset duct are 
shown in Fig. 14 for the supersonic core conditions in Table 1. The data for Ab/Ac =1.0 and Ab/A c = 2.5 are indicated 
with dashed and solid lines, respectively. For Ab/A c = 2.5, the peak frequency band for 140° at NPR t = 1.0 is 1000 Hz 
for the core and bypass exhaust conditions used in Fig. 14. The results in Fig. 14 indicate maximum noise reduction 
resulting from the addition of the offset tertiary stream occurs at frequencies above the peak frequency for the 
simulated two- stream jet. For all core and bypass conditions, maximum noise reduction occurs for NPR t = 2.1 and is 
accompanied by high-frequency noise increases. For the same bypass, core, and tertiary conditions, maximum noise 
reduction for Ab/A c = 1.0 is greater than that for Ab/A c = 2.5. Noise reduction of up to 8 dB is achieved for NPR C = 
2.3, NPR b = 1.8, and NPR t = 2.1. 



(a) 



(c) 




(d) 


Figure 14. The one-third octave spectra noise reduction on the thick side of the offset duct nozzle in the peak- 
jet-noise direction (140°) for (a) NPR C = 2.1 and NPRb = 1.8, (b) NPR C = 2.1 and NPR b = 2.1, (c) NPR C = 2.3 and 
NPR b = 1.8, and (d) NPR C = 2.3 and NPRb = 2.3. All data were acquired with NTR C = 3.0 and Mfj = 0.3. Dashed 
and solid lines are for Ab/A c = 1.0 and Ab/A c = 2.5, respectively. 


The experimental results were used to develop a one-third-octave-band noise-reduction model for acoustic 
radiation on the thick side of the jet given by 


NR\/ - Po + P \ X l +P 2 X 2 + A X 3 + Pa X A + P 5 X 2 + P 6 X 2 X 3 + P 7 X 3 X 4 +/3 8 X 2 x a 


( 1 ) 


where NR X/ is the noise reduction in one-third-octave bands and (3 ’ , n - 0, • • *, 8 are coefficients determined from the 

73 

method of least squares. The values for x l ,x 2 , x 3 and x 4 are given by the centered and standardized values of Ab/A c , 
NPR t , NPR C , and NPRb. Centered and standardized values shift and scale the quantities used in the experiments to fit 
a span of -1 to 1. For the bypass-to-core area ratio, Ab/Ac =1.0 becomes -1.0 and Ab/Ac = 2.5 becomes 1.0. The 
model, developed for the supersonic core conditions in Table 1, used centered and standardized values for the nozzle- 
pressure-ratio terms based on values for NPR t between 1.5 and 2.1, NPR C between 2.1 and 2.3, and NPRb between 1.8 
and 2.3. The last three terms in the equation represent interaction terms. The inclusion of interaction terms indicates 
that the impact of increasing the level of one variable in the interaction term on the resulting noise reduction depends 
on the level of the other variable in the interaction term. The resulting coefficients are shown in Fig. 15 for a range 
of observation angles. The values of p o indicate that noise reduction increases with increasing polar angle. Increasing 
the level of Ab/Ac (the Pi term) reduces the maximum noise reduction at large angles to the jet inlet axis (160°). 
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Increasing the level of NPR t (the P 2 term) increases peak-noise reduction at large angles to the jet inlet axis (130° - 
160°) and increases high-frequency acoustic radiation at all polar angles. Increasing NPR C (the p3 term) increases 
mid-frequency acoustic radiation in the peak direction. Increasing NPRb decreases mid-frequency radiation in the 
peak direction and increases high-frequency radiation at small observation angles. The interaction and quadratic terms 
impact noise reduction at high frequencies. 
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Figure 15. The one-third-octave-band noise reduction coefficients for the thick side of the offset duct nozzle 
system. The polar angles are indicated in each plot. 
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IV. Acoustic Analogy 


A. Basic equations 

The noise prediction method used in this paper is based on the generalized acoustic analogy of Goldstein [27], 
which has been used to develop noise prediction methods for round, cold jets by Goldstein & Leib [3], Leib & 
Goldstein [28], and for rectangular jets by Leib [29] . 

Goldstein and Leib [3] show that the acoustic spectrum at x due to a unit volume of turbulence at y is given by 


4 > (* 1 y) = ( 2 ^f r Aj (*| j;©) J Ki (*| J +J 7 (y,tj,co)dt], ( 2 ) 

V 

where the Greek indices range from one to four and the Latin indices from one to three, an asterisk denotes the complex 
conjugate, T Xj is the Fourier transform of a propagator function, and 0 % XjKl is the source spectrum. 

The propagator function is defined as 

oo 

r J e ' a{t ~ T) rAj (*| y,t-r)d(t-r), (3) 

—00 

where 


r Aj (x,t\y,r) 


Sg a M (y,T |*,f) 


(r-i)8u^8 a u{y,T\x>t), 


(4) 


with y being the ratio of specific heats, v k the Fa vre- averaged velocity, and 8m (y, r\x 9 t) the fourth component of 

the adjoint vector Green’s function, which can be calculated from equations (4.8) - (4.11) of reference [3] once the 
mean flow is specified. 

The source spectrum ^ y Jki e 4 uat i° n (2) is related to the spectrum of the generalized Reynolds stress auto- 


covariance tensor, 


with 


O V 

H ,j«i(y^v) = ^ | e~ UOT R xjKl (y,^,T)dT, 


(5) 


R xjAy ^ T ) 1 


pv '/j- pv '/j 


(y> T ») 


K v rK v 'i 


(y + Tj,T 0 +r), 


( 6 ) 


where v[ =v x - v x denotes a generalized, four-dimensional "velocity" fluctuation, with v i9 i — 1,2,3 being the 

( i \ 

ordinary fluid velocity and v A =[y- 1) h' + — v' 2 , where ti is the fluctuating enthalpy, and an overbar represents 

V 2 ) 

time average, by the simple linear transformation 




X j kI 


— p j-j p 

6 Aj ,0-m 1 1 (jmjun*' kI ,/ un ’ 


(7) 


where 


Aj,an 


1 8 A* 8 jin' 


r - 1 


s 3 s 

Aj an 


The acoustic spectrum is obtained upon integrating I m (x I y ) over the volume of the jet 


I < a ( x ) = \L( xi y) d y- 

v 
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( 8 ) 

( 9 ) 



B. Green’s function 


For the noise predictions in this paper, the mean flow is represented by a unidirectional transversely shear mean 


flow 


= S n U(y 2 ,y 3 ), c 2 =c 2 (y 2 ,y 3 ), p= constant, 


( 10 ) 


for which the adjoint vector Green’s function can be expressed in terms of a single scalar function. For observer 
locations in the far field, the relevant problem for the Green’s function can be expressed in polar coordinates as ([30]) 


d 2 g 1 dg 1 d 2 g dg 1 dg 2 
— f- + — + + ^ — + ? — + CO 2 '&g=0 , 

ctyr Jr ^3Y Jr Jr 


( 11 ) 


where 


and 


^ ( Jr ’ % ) 


?(y T ’<p 0 ) 


j 2coSi9 dM(y T ) 1 dc 2 (y T ) 

j[l-M (j r )cos <9] dy T c 2 (y T ) dy T j’ 

| 2 cos <9 SM (jv) + 1 3c 2 (j> r )) 

[l-M (j r )cos i9] y T d<p 0 c 2 (y T ) y T d % 


&(y T >w o) 


[l-M (y r )cos i9] 2 cos 2 19 
c 2 ( y T ) cl 


( 12 ) 

(13) 


(14) 


with M (y r ) = f/ (y r )/c oo being the acoustic Mach number, the ambient sound speed and c 2 (y r ) the mean 
sound speed squared, subject to the far- field boundary condition 

(cole ) 2 e ~ i ^ CooSin&y T COS i <P ~ <Po ) e i’r/4 

g(y T ,(p 0 ;(p,$:co ) -> — 2 + outgoing waves, ( 1 5) 

2(2 7l) sin i9 co / c x 

as y r = \y T \ = + y 2 —> 00 , where 3 = cos 1 (jq/x) ( x = |jc| ) is the observer polar angle measured relative to the 

downstream jet axis, cp - tan -1 (x 3 / v 2 ) the observer azimuthal angle and (p 0 - tan -1 ( J 3 / J 2 ) • Noise prediction results 

will be presented in terms of the observer polar angle measured relative to the nozzle inlet, 6 = 180 - 3 . 

For the non-axisymmetric mean flows considered in this paper, the approach of Leib [30] is used to obtain the 
Green’s function. The approach consists of expanding the coefficients in (11), and the solution, g , in terms of their 

azimuthal Fourier coefficients, approximating the transverse derivatives by second-order finite differences and solving 
the resulting banded system of algebraic equations using a sparse system algorithm. Further details of the numerical 
methods can be found in reference [30]. 


C. Source model 

The source model used in this work is the hybrid (space- time/spectral) model of Leib and Goldstein[28]. This 
model was constructed with the time and streamwise separation dependence of the Reynolds stress auto -covariance 
tensor, (6), represented by a functional form accounting for experimentally observed features of this quantity, and the 
transverse separation dependence specified in terms of a frequency-dependent length scale that enters through its 
spectrum. The approach was referred to as a hybrid one, to signify the combined used of space-time and wavenumber- 
frequency domain modeling. 

The model provides a formula for the source spectral function, 
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(16) 


' p xj*i(y’ k ’‘ a ) s S e ' 


= r '(*•*) 


as ([28]) 


XlJ\ 


lj (l + ^ 

u„ 


1/2 00 


■Z 

m,l=0 




DTD 1 

co 


i a i -k 2 T /4{i+R 2 ) m 

RdR(l + R 2 ) e 


(17) 


where a m t are constants, l t ,i = 0,1 are turbulent length scales, the frequency-dependent transverse length scale, l T , 
is modeled as 


kM 

[l + <w 2 (l + ^ 3 )/(l + Z>)] 1/4 


( 18 ) 


where b is a constant, the normalized frequency cb = col 0 /U c , with co being the radian frequency, U c the 
convection velocity of the turbulence, k 2 = ( ' k\ + k 2 ) (l + ) is a normalized transverse wavenumber and 

R 1 = or + k 2 ,with k x = {k x - co! U c )l x . 


The operators D k and D~ , are defined as 




(id d ^ 

(i A 

1 — 

— co + k. 

J 0 deb dk Xy 

l 1 

Vo y 


D n = — 

O) ~ 

CCQ 


(19) 


In [28] and [29] the source model described above was used to model the momentum- flux source terms of the 
generalized acoustic analogy formulation for noise predictions in cold jets. In these calculations, the enthalpy-flux 
source terms were neglected (ie. R Ajicl was replaced by R ijkl ). Afsar, Goldstein and Fagan [31] analyzed the enthalpy- 

flux source terms of the generalized acoustic analogy with the aim of reducing the number of independent spectral 
components contributing to the acoustic spectrum and understanding the structure of the remaining terms. For the 
former, they used the symmetry properties of the generalized Reynolds stress auto covariance tensor, introduced two 
kinematic approximations (consistent with those made in [28] and [30] for the momentum-flux source terms) and 
exploited the disparity of the turbulence correlation and mean flow length scales. Their analysis expressed the total 
acoustic spectrum in terms of contributions from momentum-flux auto covariance, enthalpy- flux/momentum- flux 
covariance and enthalpy-flux auto covariance source terms. The latter two are additional terms, beyond the 
momentum- flux terms, that generally need to be included for noise predictions in heated jets. 

In this paper, predictions for cold jets are carried out neglecting all enthalpy-flux source terms, and the 
formulation in [31] is used to include the enthalpy- flux auto covariance source terms, which are expected to be the 
dominant ones in the moderate Mach number, moderately heated jets considered here, in a case where the core stream 
is heated relative to the ambient. We use the same functional form described above (but with different source 
coefficients) to model the enthalpy- flux auto covariance source terms. An approximation for the amplitude scaling of 
these terms is obtained by assuming a quasi- normal form and neglecting enthalpy- momentum correlations as 

KiM (y,o, 0 ) - (r- if -(20) 

Since the Reynolds-averaged Navier-Stokes solver used in this work does not provide for evaluation of the enthalpy 
(or temperature) variances, an approximation for the total enthalpy variance, (h’ t j , was obtained using the Empirical 

Temperature Variance (ETV) model ([32]) . The latter was developed for single-stream round jets and its use here is 
an extrapolation of the model to multi- stream jets. 
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V. Prediction Results 


A. Flow 

In this sub-section, results from the Reynolds- averaged Navier-Stokes (RANS) solutions used as input to the noise 
prediction method are presented and compared with PIV data. Examples of the quality of the representation of the 
azimuthal variation of the mean- flow-dependent coefficients of the Green’s function equation (11) by a finite number 
of azimuthal modes in an offset third stream case are shown. 

The RANS solutions were obtained using the WIND US code ([33]) with the Mentor Shear Stress Transport (SST) 
model. 

1. Axisymmetric 

Figures 16 and 17 show contours of the acoustic Mach number and turbulent kinetic energy obtained from the 
RANS solutions for the axisymmetric case with Ab/Ac=1.0, A t /Ac=1.0, NPR c =1.8,NPRb=1.6, with all streams cold, 
with and without a third stream, respectively. Results from the RANS solutions suggest that the effect of the third 
stream is to move the end of the primary (marked by vertical bars) downstream relative to the case without a third 

stream. The levels of peak turbulent kinetic energy in the two cases are about the same, (k/u] ) ak = 0.026, with the 
peak occurring very slightly further downstream in the case with the third stream. 




Figures 18 and 19 show comparisons of results from the RANS solutions with PIV data for the case with the third 
stream. In Fig. 18(a), contours of the mean streamwise velocity, normalized by the ideal core jet exit velocity, are 
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shown, with the RANS solution in the upper part and the PIV data in the lower part. Figure 18 (b) shows the 
corresponding turbulent kinetic energy. Figure 18 (a) shows that the RANS predicts a slightly shorter primary core 
length compared with the data and Fig. 18 (b) that significantly higher turbulence levels are predicted by the RANS. 
Figure 19 shows the streamwise variation in the normalized mean streamwise velocity on the centerline and at a 
location just off the centerline, and the normalized turbulent kinetic energy at a radial location near the maximum 
turbulence level, from the RANS and PIV. On the centerline, the RANS over-predicts the length of the viscous wake 
behind the nozzle plug, while just off the centerline the RANS can be seen to predict a slightly shorter core length. 
The comparison of normalized turbulent kinetic energy shows that the RANS over-predicts the turbulence levels near 
the end of the plug, but levels closer to the data are predicted near the end of the measurement region. RANS solutions 
are known to generally over-predict the turbulence levels in high-speed jet flows ([34]). The differences here are 
somewhat greater than those seen in simpler geometries, possibly indicating additional challenges for RANS-based 
calculations in more these more complex nozzle configurations. We note that the turbulence levels obtained from the 
RANS solution for the case without a third stream (not shown) are in somewhat better agreement with the PIV data. 




V / D ® 

*1 _ equiv 


(a) 


(b) 


Figure 18. Comparisons of RANS solution with PIV data for three-stream axisymmetric case, (a) Normalized 
mean streamwise velocity, (b) Normalized turbulent kinetic energy. 



(a) (b) 

Figure 19. Comparison of RANS solution with PIV data, (a) Normalized mean streamwise velocity, (b) 

Normalized turbulent kinetic energy. 


2. Offset Third Stream 

Figure 20 shows contours of the mean acoustic Mach number and turbulent kinetic energy from the RANS solution 
in a number of axial slices through the jet, for the case: Ab/A c =1.0, A t / A c =1.0, NPR C =1.8, NPRb=1.6, NPR t =1.2 , with 
all streams cold and the third stream offset. The third- stream offset results in a thick and a thin side to the jet, in terms 
of its mean velocity profile, and a significant shift of the turbulent kinetic energy to the thin side. 

Figure 21 shows comparisons of the normalized mean streamwise velocity from the RANS solutions and PIV data 
for this cold offset case. The plots show contours of the normalized mean streamwise velocity in cross-flow planes at 
a number of streamwise locations downstream of the end of the nozzle plug. Experimental data is shown in the upper 
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parts and RANS results in the lower parts of these plots. Close to the end of the nozzle plug there are significant 
differences between the shapes of the contours from the PIV and RANS. Further downstream the shapes become more 
similar, but the RANS solution seems to mix out to an axisymmetric mean flow closer to the end of the nozzle plug 
than the data indicates. Comparisons (not shown) of turbulent kinetic energy from the RANS solution with PIV data 
for this case show that the former predicts the significant shift of higher turbulence levels to the thin side of the jet 
as seen in the data (see Fig. 20 (b)), but the peak levels exhibit differences similar to those of the axisymmetric case 
in the last sub- section. 



(a) (b) 

Figure 20. Contours of (a) mean acoustic Mach number and (b) turbulent kinetic energy. 
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Figure 21. Comparison of the normalized mean streamwise velocity in cross-flow planes at a number of axial 

locations. Cold, offset third stream case. 
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Figure 22 shows comparisons of the normalized mean streamwise vorticity from the RANS solutions and PIV 
measurements. The presentation of these plots is the same as those in Fig. 21, but different streamwise locations are 
shown. Very close to the end of the nozzle plug, a region of significant mean streamwise vorticity was found in the 
PIV data. The RANS solution also exhibits this streamwise vorticity, but at a slightly different location and the RANS 
result tends to dissipate the vorticity too quickly as the flow evolves downstream. The generation of significant levels 
of streamwise vorticity in the region near the end of the nozzle plug would result in enhanced mixing and could be 
partially responsible for the reduction of low- and mid-frequency noise observed in some of the three-stream nozzle 
configurations studied. The model used for noise predictions in this paper does not explicitly account for the presence 
of mean streamwise vorticity (see equation (10)), but some of its effects, through its impact on the mean streamwise 
velocity and turbulence, may be manifest indirectly. 



Vi! Deauiv — 1 yi/D eauiv 


Figure 22. Comparisons of normalized mean streamwise vorticity. 


Figure 23 shows the azimuthal variation of the mean-flow-dependent coefficients of equation (11), fori9 = 30 ( 
0 = 150 ) , determined from the RANS solution and as represented by four azimuthal Fourier modes at two locations 
within the jet. There is some discrepancy between the Fourier representation of R near cp Q - 0, n , but for the most 
part, the coefficients are well represented by this relatively small number of modes, an important result for the viability 
of the noise prediction method for non-axisymmetric jets. 
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Figure 23. Azimuthal variation of the mean-flow-dependent coefficients of equation (11) and their 
representation by four Fourier modes. R , red; ? , green; ^ , blue. Curves from RANS solution, symbols 
from Fourier series, (a) yjD equiv = 2. 0, y T /D equtv = 0. 07; (b) y 1 /D equiv = 4. 0,y T /D equlv = 0. 4. 

B. Noise predictions 

1. Axisymmetric 

In this subsection, prediction results are presented, based on the method in reference [28], for two axisymmetric 
cases: Ab/Ac=1.0, A t /A c =1.0, NPR c =1.8,NPRb=1.6, with all streams cold, with (NPR t = 1.2) and without (NPR t = 1.0) 
a third stream. The coefficients in the source model are the same as those used in [28]. 

Specifically, for the amplitudes of the various components of the Reynolds stress auto-covariance, 
R ijkl (y,0,0) , it is assumed that 

R ijU (y,0,0)*C ijkl (pk) 2 , (21) 

where the C jjkl are constants, p denotes the mean density and k the turbulent kinetic energy determined from the 
Reynolds- averaged Navier-Stokes (RANS) computation. As in [28], C 2222 = C 3333 = 0.159C 11U 

C 2 233 = 0.0, C n22 = C 1133 = 0.047 C nn and C 1212 = 0.288C im . In order to set the absolute level of the turbulence (and 
therefore of the sound), the value for axial component of (21) is set as C im - 2^v[ 2 j / k 2 where, v' = v 2 , v'} , and 

v' 2 =0.8 k . The length scales, ,/ 2 = / 3 =l T , are related to the RANS computations by setting l x = C x ^ / and 

jV 2 / 

l T —C T K / e • The specific values used for the constants in these length scales are shown in Table 2. The values of 
a n o / a o o are taken to be independent of the source location y and the values used for the ratios of the coefficients 

in the truncated series representation (17) are indicated in Table 3. In the model (18) for the transverse frequency- 
dependent length scale, b = 0.1 for the quadrupole-like terms and b - 0.5 in the dipole-like term (see [3], [28]). As in 
[28], U c - 0.68f/ (/ where U d is the jet centerline velocity, for the quadrupole-like source terms and U c - 0.8 U cl in 
the dipole- like term. 

Table 2: Constants used to determine the length scales for the indicated components of the turbulence auto- 
covariance 


Component 

Co 

C, 

C T 

1111 

0.7 

1.2 

0.4 

2222=3333 

0.7 

0.8 

0.89 

1122=1133 

0.7 

1.1 

1.0 

1212 

1.05 

1.0 

1.1 
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Table 3: Coefficients in the truncated series representation of the turbulence spectrum. 


Component 

a i,o / 
/ a o,o 

a 2,0 / 
/ a o,o 

a 3,0/ 

/ a o,o 

1111 

0.073 

0.070 

-8.48x10-4 

2222=3333 

0.519 

0.049 

-0.0097 

1122=1133 

0.103 

0.079 

0.0 

1212 

0.559 

-0.006 

-0.015 






Figure 24 shows comparisons of the prediction results using 
this model with experimental data for the cold, axisymmetric case 
with NPR t =1.2. Reasonably good agreement is obtained, with the 
predictions being generally within about 2 dB of the data. 

Figure 25 shows the effect of the third stream on the noise for 
the cold, axisymmetric case. The experimental data (shown offset 
by 10 dB relative to the predictions) shows that the introduction of 
the third stream at NPR t =1.2. results in a relatively small (about 
1 dB or so) increase in noise at low-frequencies at all polar angles 
shown, and a slight noise reduction at very high frequencies at 
ninety degrees. The calculations predict that the introduction of the 
third stream increases the noise levels by a little more than 1 dB 
across nearly all frequencies at ninety degrees, with less effect at 
high-frequencies in the downstream polar angles. 

Figure 24. Noise predictions for a cold, 
axisymmetric three-stream jet. Curves 
offset by 5 dB for clarity. 




Figure 25. Effect of third stream on axisymmetric cold case. 


2. Offset Third Stream 

In this subsection, comparisons of prediction results are presented, using the approach in reference [30] to obtain 
the Green’s function, with experimental data for the case: Ab/A c =1.0, A t / A c =1.0, NPR C =1.8, NPRb=1.6, NPR t =1.2 , 
with all streams cold and the third stream offset using an offset duct. The specific values used for the coefficients of 
the RANS-based length scales and of the ratios of the coefficients in the truncated series representation (17) are the 
same as those used for the non-circular cases in [29]. The former are in listed in Table 4 and the latter are the same 
as in the axisymmetric calculations (Table 3). In the offset jet predictions, the convection velocity, U c , is taken to be 
equal to the local mean velocity at the point of maximum turbulent kinetic energy. 
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Table 4: Constants used to determine the length scales for the indicated components of the turbulence auto- 

covariance 


Component 

Co 

C, 

C T 

1111 

0.7 

0.95 

0.32 

2222=3333 

0.7 

0.63 

0.71 

1122=1133 

0.7 

0.87 

0.79 

1212 

1.05 

0.79 

0.89 



Figure 26 shows the results of the comparisons for the cold 
offset case. The peak noise level at 90 degrees, and its variation 
with polar angle in the downstream direction is relatively well 
predicted, as is the shape of the 90-degree spectrum. At 150 
degrees, the predicted results peak at a lower frequency than the 
acoustic data, and do not exhibit the same degree of azimuthal 
variation between the thick and thin sides of the jet. It is possible 
that the rapid mixing of the flow in the RANS solution is partially 
responsible for the relative lack of azimuthal directivity in the 
noise predictions 



3. Hot jet predictions 

In this subsection comparisons of prediction results with experimental data are presented for the case: Ab/A c =1.0, A t / 
A c =1.0, NPR C =1.8, NPRb=1.8, NPR t =1.2 , NTR C = 3.0, NTRb = NTR t = 1.25, with the third stream offset. For these 
calculations, the coefficients in the momentum-flux source terms were the same as those used in the cold offset 
predictions of the previous sub-section. The coefficients in the enthalpy-flux source terms (41,41 and 42,42) were 
taken to be the same as those of the corresponding momentum- flux terms (11,11 and 22,22), except for the coefficients 
of the transverse length scales, C T . The latter were set to: C T - 1.5 for the 41,41 component and C T =0.8 for the 
42,42 component. It was found necessary to increase the values of these coefficients above those used in the 
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corresponding momentum- flux terms to match the rapid increase of the peak noise levels with increasing polar angle 
in the heated case. The scaling of the amplitude of the enthalpy-flux source terms was set using an approximation 
based on the Empirical Temperature Variance Model of [32]. 

Figure 27 shows comparisons of the predictions with data for the heated case. Included on these plots are results 
obtained using only contributions from the momentum-flux source terms. These results are significantly below the 
data at all polar angles, suggesting the need for the additional enthalpy-flux source terms in this case. The results 
including the enthalpy- flux source terms are much closer to the data, and match the rapid increase of the peak noise 
levels with increasing polar angle in the heated case reasonably well. As mentioned above, the coefficients of the 
transverse length scales in the enthalpy-flux source terms were increased, making these terms more directional, to 
obtain these results. 



Figure 27 - Comparisons of noise predictions with data, hot offset case. 


VI. Conclusions and Discussion 

The introduction of an axisymmetric tertiary stream to a dual-stream nozzle system reduces acoustic radiation in 
the peak-jet-noise direction. The largest noise reductions (up to 3 dB) are achieved for a nozzle system with a small 
bypass-to-core area ratio (Ab/Ac =1). For small and broadside angles to the jet, the introduction of the third stream 
has no impact on low- and mid-frequency radiation for a large bypass-to-core area ratio (Ab/A c = 2.5) or for a small 
tertiary-to-core area ratio (A t /Ac = 0.6) and increased low-frequency radiation for a nozzle system with Ab/A c =1.0 
and A t / A c = 1.0. An inverted third-stream velocity ratio, where the velocity of the third stream is larger than that of 
the bypass stream, results in elevated high-frequency radiation for all nozzle area-ratios and observation angles. For 
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the same operating conditions and nozzle configuration, greater noise reduction is achieved for static simulated 
forward- flight conditions than for Mfj = 0.3. 

The introduction of asymmetry to the tertiary stream through partial blockage of the tertiary nozzle exit using nozzle 
inserts provides no acoustic benefit. Introduction of asymmetry to the tertiary stream using an offset duct placed just 
upstream of the tertiary nozzle entrance results in noise reduction (relative to the simulated two- stream jet) in the peak 
jet-noise direction on the thick side of the jet. For small bypass-to-core area ratios (Ab/Ac = 1.0), the introduction of 
the third stream at NPR t >1.3 results in increased acoustic radiation throughout the spectra at small and broadside 
polar angles on the thick side of the jet and at all polar angles on the thin side of the jet. The offset duct is more 
effective at noise reduction for supersonic core conditions than for subsonic core conditions with a maximum noise 
reduction of 8 dB in the peak-jet-noise direction. At supersonic core conditions, noise reduction in the peak-jet-noise 
direction increases with increasing core nozzle-pressure ratio and decreasing bypass-to-core area ratio. 

A small subset of the experimental cases was chosen for use as test cases for a RANS-based noise prediction 
scheme. RANS solutions were obtained for these cases and results were compared with PIV data. The RANS solutions 
were generally able to predict the trends associated with the effects of the third stream on the flow, but quantitative 
differences were found. It is possible that, for these relatively complex nozzle systems, the use of more advanced 
turbulence models could result in better predictions of the flow. Within the RANS context, an algebraic Reynolds 
stress model has been shown to provide improved flow predictions, compared with a two -equation model, for a high- 
speed elliptical jet with relatively little additional cost ([35]). An algebraic Reynolds stress model would also provide 
additional information about the flow, such as the degree of anisotropy of the turbulence, which could be used to 
inform source models for noise predictions. Alternatively, a Large Eddy Simulation (LES) of the flow field could be 
carried out. The latter, of course, involves significantly higher computational costs, but the resulting unsteady flow 
solutions could provide, in addition to improved single-point turbulence statistics, further features of the flow, such 
as two-point, time-delayed velocity and velocity-enthalpy correlations and associated convection velocities. The use 
of more advanced turbulence models and LES for flow and noise predictions in three-stream jets will be a topic for 
follow-on research. 

Noise predictions were calculated for the selected test cases using an acoustic analogy-based formulation, with the 
RANS solutions as input. The prediction method combines numerical solutions for the Green’s function of the residual 
equations of the acoustic analogy in the locally parallel mean flow approximation with a source model constructed 
partially in the time-space domain and partially in the frequency- wavenumber domain. Results for an axisymmetric 
three- stream configuration, with all streams unheated, were in reasonably good agreement with data, but tended to 
over-predict the effect of the third stream on the noise at high frequencies for polar angles near ninety degrees to the 
jet axis. 

Results for an offset third stream, also with all streams unheated, tracked the roughly 12 dB increase in noise from 
polar angles between ninety degrees and the peak noise direction, but over-predicted the shift of this peak to lower 
frequencies. Also, although the model qualitatively predicts the experimentally observed lower noise for observer 
locations on the thick side of the jet relative to the thin side at polar angles near the peak noise direction, the amount 
of azimuthal directivity is under-predicted. Reasons for this may have to do with the relatively quick mixing and 
transition to a nearly axisymmetric mean flow in the RANS solutions, or inadequacies of the source model in 
representing the azimuthal variation of the turbulence. Work is currently under way on an acoustic analogy 
formulation, and accompanying source modelling, that treats the azimuthal coupling of propagation and source 
components in a more rigorous way. 

Results for an offset third stream with the core stream heated relative to the ambient showed the importance of 
including enthalpy- flux source contributions in noise predictions for heated jets. Predictions for this case using only 
the momentum-flux source terms of the generalized acoustic analogy were more than 5 dB below the data for polar 
angles in the downstream quadrant. Inclusion of contributions from the enthalpy-flux auto covariance source terms 
gave much improved agreement with the data. It was found that model coefficients for the enthalpy- flux auto 
covariance source terms needed to be chosen to make the contribution from these terms more directional than the 
corresponding momentum- flux terms in order to match the rapid rise in noise level with increasing polar angle for the 
heated case. Work is also under way on improving the modeling of the enthalpy-flux auto covariance source terms 
and incorporating the momentum- flux/enthalpy- flux covariance (coupling) terms in heated jet predictions. 
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